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ABSTRACT 

Prior to the explosion of a carbon-oxygen white dwarf in a Type la supernova there is a long "simmering," 
during which the '^C + '^C reaction gradually heats the white dwarf on a long (~ 10^ yr) timescale. Piro & 
Bildsten showed that weak reactions during this simmering set a maximum electron abundance at the time 
of the explosion. We investigate the nuclear reactions during this simmering with a series of self-heating, at 
constant pressure, reaction network calculations. Unlike in AGB stars, proton captures onto ^-Ne and heavier 
trace nuclei do not play a significant role. The '-C abundance is sufficiently high that the neutrons preferentially 
capture onto '-C, rather than iron group nuclei. As an aid to hydrodynamical simulations of the simmering 
phase, we present fits to the rates of heating, electron capture, change in mean atomic mass, and consumption 
of '^C in terms of the screened thermally averaged cross section for '~C + '^C. Our evaluation of the net 
heating rate includes contributions from electron captures into the 3.68 MeV excited state of '^C. This results 
in a slightly larger energy release, per '^C consumed, than that found by Piro & Bildsten, but less than that 
released for a burn to only -"Ne and -^Na. We compare our one-zone results to more accurate integrations 
over the white dwarf structure to estimate the amount of '^C that must be consumed to raise the white dwarf 
temperature, and hence to determine the net reduction of during simmering. 

Subject headings: nuclear reactions, nucleosynthesis, abundances — supernovae: general — white dwarfs — 
galaxies: evolution 



1. EMTRODUCTION 

The leading scenario for a Type la supernova (hereafter 
SNe la) is the thermonuclear incineration of a carbon-oxygen 
white dwarf that has increased in mass, through accretion, 
to just below the Chandrasek har limit (for a review, see 
iHillebrandt & Niemeyej |2000() . Despite the advances in 
modeling the post-i gnition flame evolution (for a sampling 
of recent work, see iGamezo etal.ll2004t iPlewa etalj [20041: 
iRopke et alj2006tlJordan et al.l2007hT we still lack a complete 
understanding of which subset of the binary white dwarf pop- 
ulation become SNe la, and how differences in the progenitor 
map onto differences in the outcome of the explosion. 

The composition of the white dwarf at the time of the 
explosion should have an effect on the nucleosynthesis that 
takes place during the explosion and the isotopic abundances 
of the final composition. The observable properties of 
SNe la resulting from Chandrasekhar-mass explosions are 
chiefly determined by their final composition, the velocity 
profiles of key spectral lines at early- and late-times (e.g., 
P-Cygni profiles of Si II at 615.0 nm at early times), the 
opacity of the material through which the photons from 
radioactive decay must propagate, the kinetic energy of 
ejecta, and its interaction with the density profile of the 
surro un ding circumstellar o r inter st ellar medium ( Filippenl^ 
119971: iPinto & Eastmaiil 120001: iHiUebrandt & NiemeveJ 



20001 iLeibundsutI 120011 iMazzah & Podsiadlowskil 


2006 


Marion et alj 20061 


Blondin et al.ll2006h iBadenes et al.1 


2007 


Wooslev et al.l 12007 


). The dominant parameter in setting the 



Electronic address: chamulak@pa.ms u.edu| 
Electronic address: ebrown@pa.msu.e3u] 
Electronic address: fxt44@mac.com 
Electronic address: dupczakk@msu.edu 

' Department of Physics and Astronomy, Michigan State University, East 
Lansing, MI 48824 

^ Joint Institute for Nuclear Astrophysics 

' National Superconducting Cyclotron Laboratory 

* Thermonuclear Applications, X-2, Los Alamos National Laboratory 

' School of Earth and Space Exporation, Arizona State University, Tempe, 
AZ 85287 



peak brightness, and hence width, of the light curve is widely 
believed to b e the m ass of ^*Ni ejected by the explosion. 
iTimmes et al.l (l2003h showed the mass of ^^Ni produced 
depends linearly on the electron fraction, Yg, at the time 
of the explosion, and that Y^, itself depends linearly on the 
abundance of ^-Ne in the white dwarf. 

In this paper we explore, using a reaction network coupled 
to an equation for self-heating at constant pressure, the re- 
duction in Ye that occurs after the onset of the thermonuclear 
runaway (when the heating from '-C + '~C reactions is faster 
than cooling by thermal neutrino emission), but before the 
burning becomes so fast that local regions can thermally run 
away and launch a flame. This "simmering" epoch lasts for 
~ 10^^ yr, long enough that electron captures onto products 
of '^C burning can reduce the free electron abundance Ye- A 
similar, but independent calculation , of the reduction of Y e 
during simmering was performed bv lPiro & Bildsteiil (l2008l) . 
Their calculation did not use a full reaction network, but it 
did take into account the change in energy of the white dwarf 
due to the growth of the convective zone. Our calculation 
agrees with their findings, in particular that there is a maxi- 
mum Ye at the time of the explosion, and that the reduction 
in Ye is linear in the amount of '~C consumed prior to when 
the rate of '^C burning outpaces that of the weak reactions. 
This paper expands on their work in three ways. First, by us- 
ing a full reaction network, we are able to quantify the role of 
-^Ne and trace nuclides in setting the change in electron abun- 
dance with '^C consumption, dYeldYn- Second, we calculate 
the heating from the electron capture reactions and include 
the contribution from an excited state of '^C. Third, we pro- 
vide tabulated expressions for the rate of heating e, the rate of 
change in electron abundance dYe/dt, and the rate of change 
in the mean atomic mass d{A}/dt in terms of the reaction rate 
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for the '^C + '^C reaction. These expressions are useful input 
for large-scale hydrodynamical simulations of the simmering 
phase which do not resolve such microphysics. 

We first give, in § |2] a simple estimate for the reduction 
in Yt, during the pre-explosion simmering and describe the 
role of ^^Ne and other trace nuclides. In §[3] we describe our 
numerical calculations, explain the reactive flows that occur 
(§ 13. Il l, and give simple approximations to the heating rate 
and carbon consumption (§ 13. 2b . We detail, in § 13.31 some of 
the limitations of our approach. We evaluate the energy re- 
quired to raise the white dwarf central temperature, and hence 
the amount of '^C that must be consumed, and compare it 
against the one-zone calculation. We also estimate the cen- 
tral temperature at which convective mixing becomes faster 
than electron captures. This convective mixing advects elec- 
tron capture products to lower densities where they can /3^- 
decay: the convective Urea process. Because each electron 
capture-decay cycle emits a neutrino-anti-neutrino pair, there 
is energy lost from the white dwarf, and our calculation under- 
estimates the amount of consumed prior to the flame run- 
away. The convective Urea process (Paczyhski 1972]) reduces 
the rate of heating by nuclear reactions (thereby increasing 
the amount of '^C that must be consumed to raise the temper- 
ature), but cannot result in a net decrease in entropy and tem- 
perature for constant or increasing density jStein et al.ll 19991: 
[Stein & Wheeler 2006). The Urea reactions also tend to re- 
duce the effects of buoyancy, and in degenera te matter have 
a dire ct influence on the convective velocity (iLesaffre et al.l 
l2005h . The paper concludes (§ |4|i with a discussion of the 
implications for systematic variations in the mass of ^^Ni pro- 
duced in the explosion. 

2. THE REDUCTION IN ELECTRON ABUNDANCE DURING THE 
EXPLOSION 

The demise of an accreting white dwarf begins when the 
central temperature and density are such that the heating from 
the '^C + '^C reaction becomes greater than the cooling from 
thermal neutrino emission. For a density p - 2. x 10^ gcm'-^ 
this re quires a temperature T k 3.0x 10 K (see lGasques et al.l 
I2005L for a recent calculation). Initially the heating timescale 
is long, tu = T{dT Idty^ ~ 10^ yr; as the ter nperature rises 
and th e reaction rate increases, tu decreases. IWooslev et al.l 
(120041) estimate that when T > 7.6 x 10^ K, fluctuations in the 
temperature are suflicient to ensure that a local patch can run 
away and the flame ignites. 

The basic reactions during '^C burning were first worked 
out in the context of core carbon burnin g in evolved stars 
(iReeves & SalpeterllI959l: ICamerorilll959l) . During simmer- 
ing, '^C is primarily consumed via '^C('^C, a)^''Ne and 
'^C('^C,/7)^"^Na. These reactions occur with a branching ra- 
tio 0.56/0.44 for T < lO'' K (jCaughlan & Fowler 1988). 
At temperatures below a; 7 x 10^ K, neutronization — 
that is, a reduction in Ye — occurs via the reaction chain 
^^C{p,jy^N(e^,Vt,y^C. This electron capture implies that 
there is a max imum Y^ at the time of the explosion, as first 
pointed out bv iPiro & Bild sten (2008). One can readily esti- 
mate the change in electron abundance, AFg. For simplicity, 
take the branching ratio for '^C + '^C to be 1:1 for producing 
p + ^"^Na and ''He + ^°Ne. Thus there is one p produced for 
every four '^C consumed. Two additional '^C are consumed 
via '^C(«, y)'-'CC'He, «)'^0, which also destroys one '*He nu- 
cleus, and via '^C(/7, 7)'-'N(e", Vf)'''C, which also destroys 
one p. Thus for every 6 '^C consumed there is one electron 
capture, so that dYddYn ~ 1/6, where Yn is the molar abun- 



dance of '^C. As an estimate for the heating from this reaction 
sequence, summing over the g'-values for the strong reactions 
gives a net heat release of ( 1 6 Me V) 16-2.1 MeV per '^C nu- 
cleus consumed. Note that at densities above 1 .7 x 10^ g cm"-', 
the reaction ^''Na(e", Vf)^''Ne contributes to the rate of de- 
crease in Ye, so that dYddYn ~ 1/3 at those densities. The to- 
tal effective rate of dYeldY\2 depends on the rate of convective 
mixing and the size of the convective core (see § 13. 3b but is al- 
ways at least as large as the contribution from '■'N(e", Ve)'"'C. 
In the following sections, we investigate these reactions in de- 
tail. 

2.1. The role of neon-22 and other trace nuclides 

Reactions on ^^Ne, ^^Na, and other trace nuclides also occur 
during shell-burning in asymptotic giant branch (AGB) stars, 
and we briefly summarize their role in that context before de- 
scribing the very different environment in a simmering white 
dwarf core. In AGB stars more massive than about 4 Mq, the 
hydrogen burning shell, at a temperature of 60-100 MK, ex- 
tends into the convective envelope. The envelope composition 
is then directly affected by the various hydrogen-burning cy- 
cles: CNO, N eNa, and MgAl (iLattanzio & Boothrovdl[l997l: 
lHerwigll2005l) . The hydrogen-burning shell is also disturbed 
by thermal pulses due to ignition of the helium layer. At each 
pulse, dredge-up of material may occur, in which helium- 
burnt material is mixed into the stellar envelope, pollut- 
ing it with "^He, | ^C, ^^Ne, and heavy i-process elements 
dlzzard et al.ll2007b . Thus, the fate of Ne is either to con- 
tribute to hydrogen burning via ^^Ne{p, 7)^''Na or to become 
a neutron source for the i-process via the ^^Ne{a, «)^^Mg re- 
action. The ^^Ne(Q;, n)^^Mg reaction requires the high tem- 
peratures (r > 2.5 X 10^ K) that can be found at the bot- 
tom of the pulse-driven convective zone during the helium 
shell flashes. The neutrons are release d with high density 
[log,(N„/cm-2)_~ 10] in a short burst (iGallino et al.l 119981: 
Busso et al. 19991) . These peak neutron densities are realized 
for only about a year, followed by a neutron density tail that 
lasts a few years, depending on the stellar model assumptions. 
These neutrons are the genesis of the classic high-temperature 
i-process in AGB stars. 

Should ^^Na be present, there are two usual possibili- 
ties for subsequent nuclear processing in AGB stars via ei- 
ther the ^^Na(p,a)^'^Ne reaction or the ^^Na(p,y)^'^Mg re- 
action. Th e former reaction gives rise to the classical 
NeNa cycle ("Marion & Fowledll957t iRolfs & Rodnevlll988l: 
[Rowland et al. 2004), whereas the competing (p, y) reaction 
transforms ^''Na to heavier isotopes and bypasses the NeNa 
cycle. How much material is processed through the {p,a) 
reaction on ^^Na as opposed to the competing {p,j) reac- 
tion is of interest for AGB star (and classical novae) nucle- 
osyn thesis. New measure ments of the ^^'Na + p cross sec- 
tion (iRowland et al.ll2004 suggest that for T = (20-40) MK 
^^Na(/7,7)^''Mg competes with the {p, a) branch, disrupts the 
NeNa cycle, and produces a flow into the MgAl hydrogen 
burning cycle. 

Caution about intuition developed for reaction sequences 
on ^^Ne and ^^'Na in AGB star environments seems prudent 
when applied to the dense, carbon-rich environments of white 
dwarfs near the Chandrasekhar mass limit. During the slow 
simmering preceding the explosion, the large Y12 en- 
sures that p liberated by the '^C('^C,p)^^Na branch will cap- 
ture preferentially onto '^C, rather than ^^Ne or ^^Na. Fig- 
ure [1] shows the ratio of thermally averaged cross-sections, 
A 5 A^(crv), to that for '^C(;?,y)'^N for three reactions: 
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Ne(/7, y)^-^Na {solid line), "Naip, a)^"Ne {dashed line), and 
^^Ne(/7, «)^^Na (dotted line). In addition to the Coulomb pen- 
etration, there are numerous resonances that determine how 
the cross-sections change with temperature. When the ratio 
of the thermally averaged cross-sections is of order unity, as 
it is for "Neip, y)^^Na, then the proton capture is determined 
by the relative abundances of ^^Ne and '^C. 

Note that for the latter two reactions, we plot the largest 
proton-consuming branch rather t han the (p, y) branch . We 
include screening in all reactions dYakovlev et alj|2006l) with 
the plasma taken to consist of '^C and '^O with mass fraction 
0.3 and 0.7, respectively. For T > 1 0^ K, ^[^^Ne(p,y )-^Na] 
is well-constrained experimentally (llliadis et al.ll200ll) . For 
a white dwarf with a central density p = 2 x 10^ gem"-* 
the ignition temperature (where nuclear heating dominates 
over cooling via ther mal neutrino emission) is * 3 x 10** K 
jGasques et al.l l2005h : the burning timescale becomes less 
than the timescale for -^Na(e", Vf)-^Ne once the tempera- 
ture increases beyond T > 6 x 10^ K. Over this range, 
the thermally averaged cross-sections for ^^Ne(p,y)-^Na and 
'^C(/9, y)'-'N are comparable, but the abundance of is 
far greater; having p capture preferentially onto ^"Ne would 
therefore require it to have a mass fraction roughly twice that 
of i^C. 
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Fig. 1. — Ratio of thermally averaged cross-sections, A = A^(o"v), to that 
for '^C(p, y), for thi'ee reactions: ^^Ne(p, y) {solid line); ~^Na{p, a) {dashed 
line); and ^'Ne(p, n) {dotted line). Screening is included in A; we evaluate the 
ratio at p = 3 X lO' g cm"^ for ^^Ne(p, n) and at lO' g cm"^ for the other two. 
In the temperature range where the heating timescale is slow enough for weak 
interactions to reduce (cf. Fig.|4), the thermally averaged cross-sections for 
^-Ne{p, y) and '~C(p, y) are of similar magnitude. For --Ne to compete with 
'-C for p-captures at T > 3 X 10* K requires a ^^Ne mass fraction X22 5 
(22/12)Xi2. At densities less than the electron capture threshold on ^^Na a 
small flow of ''Na(p, (}-)^''Ne can occur (cf. Fig. [2). At higher densities ^^^Na 
electron captures to form ^'Ne; the large cross-section for ^^^Ne(p, n)^^^Na 
allows it to compete with captures onto '^C if Yj^jYii > \%. 



J cm , the reaction 
•^Ne nucleus for ev- 



At densities greater then 1.7 x lO** 
^■'Na(e", V(.)^''Ne produces roughly one 
ery six '^C nuclei consumed. The screened thermally aver- 
aged cross-section for '^^Ne,{p, n)^-'Na is > 30 times that of 
'^C(p,y)'"^N (Fig. [Til, so that ^^Ne could become a compet- 



itive sink for protons. For our self-heating burn (see § 13. Il l 
starting at p = 3 x 10^ g cm""*, the abundance of ^^'Ne reaches 
723 = 3 X 10""^ (^23 ~ 0.015}'i2) by the point the heating 
timescale Th becomes shorter than the timescale for elec- 
tron captures onto ^-^Na. Although our one-zone approxima- 
tion overestimates the amount of that must be consumed 
to raise the central temperature (see § I3.3l l, should enough 
of the convective core lie above the threshold for electron 
capture onto ^■'Na, the abundance of ^■'Ne can become large 
enough to choke off t he '^C(p,y)'^N reaction, as noted by 
iPiro & Bildstenl (|2008|) . 

2.2. The reaction '^N(e", Ve)'^C 

As described above, the large '^C abundance ensures that 
protons produced by '^C('^C,/5)^^Na lead to the formation 
of /3+-unstable '^'N via '^C(/?,y)'''N unless an appreciable 
abundance of ^^Na or ^^Ne can build up. At these densi- 
ties, the rate for electron capture is substantially greater than 
the rate of yS^-decay for '^N. The electron Fermi energy is 
K 5.1 MeV(py,/10'^ gcm-3)l/^ and the ^- value for the yS+ de- 
cay of '^N is 2.2 MeV. As a result, there are several excited 
states of '^'C into which the electron can capture. Of these, 
the transition to the excited state Eexc = 3.68 MeV, with spin 
and parity J" - 3/2", is an allowed Gamow-Teller transition. 
We computed the electron capture rate usi ng the experimen- 
tal lo g /f for the ground-state transition (lAizenberg-Selovg 
Il991h . Gamow-Teller strengths to excite d states were calcu- 
lated with the shell-model code OXBASH jBrown et al.l2004l) 
emplo ying the Cohen-Kurath II potential dCohen & KurathI 
1967) in the p-she ll model space. A quenching factor of 0.67 
(Chou et aUdU) was applied to this strength, and the re- 
sulting ft values were used with the analytica l phase space 
approximations of iBecerril Reves et all (l2006h to obtain the 
capture rate. These shell-mod el calculations agre e well with 
recent (^He, f) scattering data (IZe gers et al.ll2007l) . At pY^ - 
10^ g cm""*, captures into the excited state at iSexc - 3.68 MeV 
account for > 0.3 of the total rate (/?ec = 12 s"'); this 
fraction increases with density. Because R^c > (Gp) 
this capture does not freeze out during the simmering, un- 
like capture onto ^^'Na. Although the capture into the ex- 
cited level does not increase the capture rate substantially be- 
yond that for the ground-state-to-ground-state transition, it 
does increase the heat deposited into the white dwarf from 
this reaction. We find the heat deposited from this reaction, at 
pY, = 10'' gcm-\ to be 1.3 MeV. 

2.3. Production and subsequent captures of neutrons 

Finally, we consider the contribution from heavier nuclei, 
such as ^*Fe, inherited from the main-sequence star. In 
AGB stars the reaction '''C(q', n)'^0 (during a He shell flash, 
^^Neia, «)^^Mg also contributes) acts as a neutron source for 
the i-process. In contrast, the large '^C abundance of the 
white dwarf core prevents a strong i-process flow during the 
pre-explosion simmering. The cross-section for '^C(«,y)'^C 
is 63.5 times smaller than the cross-section for ^^Fe(«, y)^^Fe 
dBao & Kappeler 1987) at an energy of 30 MeV, which is not 
sufficient to overcome the vastly larger abundance of '^C nu- 
clei (for a progenitor with solar metallicity, the '^C:^^Fe ratio 
[for a i^C mass fraction of 0.3] is 1400:1). 

3. REACTION NETWORK CALCULATIONS 

In this section we investigate the reactions that occur during 
simmering in more detail using a "self-heating" reaction net- 
work. Under isobaric conditions the temperature T evolves 
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with time according to 



dT 
dt 



s 



(1) 



Here Cp is the specific heat at constant pressure, and the heat- 
ing rate s is given by 



dt 



(2) 



where M, and F, are, respectively, the atomic mass and molar 
abundance of species /, fi^ is the electron chemical potential, 
Sy is t he neutrino loss rate, per unit mass, from the weak reac- 
tions dFuUer et al.ll982l:lLanganke & Martmez-Pin edo 200 1}, 
and A^A = 6.022 x 10^^ g"' = (1 amu)"'. We neglect 
thermal neutrino emission processes; this is an excellent ap- 
proximation over most of the integration. Our reaction net- 
work incorporated 430 nuclides up to ^^Ge and is the same 
one used t o explore the effect of Ne on the laminar flame 
speed (see IChamulak et al.ll2007l and references therein, for 
a description of the microphysics). At conditions of p = 
(1-3) X lO** gem-"* and r = 5 X 10** K, the specific heat Cp 
is dominated by the ions, which are in a liquid state (plasma 
parameter r = {Z^l'^){e^ lkYiT){pY,Np,f'^ ~ 10), and have 
Cp ~ ISk^N^I {A), where {A) is the average atomic mass. 

3.1. The reactive flows 

In this section we refine our estimate of dY^ldYn made in 
§ 12.21 We integrate equations ([T]i-(|2]i starting from the tem- 
perature at which heating from the '^C + '^C reaction equals 
the heat loss from thermal neutrino losses (this determines the 
onset of thermal instability). For simplicity, we split the solu- 
tion of the thermal and network equations. That is, for each 
time-step dt we solve the thermal equations to obtain T and p, 
integrate the reaction network at that T and p to compute Yj 
and e, and use e to advance the solution of equation ([T]i. 

In the initial phases of the simmering, the convective 
timescale is slow, and our one-zone calculations give an ade- 
quate description of the heating (when corrected for the gra- 
dient in temperature). As the temperature of the white dwarf 
increases, the heating timescale fe decreases; moreover, the 
convective mixing becomes more rapid, and one must treat 
the hydrodynamical flows in order to calculate the nucleosyn- 
thesis properly (see § 13. 3t . In this section, we will restrict 
our integration to where T < 0.6 GK, for which the heating 
timescale is > 10"^ s, so that electron captures onto ^^'Na are 
not frozen out. 

To explore the reaction channels that link '^C consumption 
with the reduction in Y^, we define the reactive flow between 
nuclides ; and j as 



dYi 
dt 



dt. 



(3) 



where the integral is for the reactions linking nuclide ; with j. 
We differentiate between F{i j) and F{j — > /), i.e., we treat 
inverse rates separately. Figures |2] and |3] show the reactive 
flows forp = 10'' gcm"^ and 3 x 10'' gcm"^, respectively. In 
both cases the composition is {X\2,X\^,X22] - (0.3,0.7,0.0). 
Each row of the chart is a different element (Z), with the 
columns corresponding to neutron number For viewing ease, 
we only plot those flows having F > 0.01 • max{F), and we 
indicate the strength of the flow by the line thickness. We 
highlight the flows of the initial '^C fusion reactions with a 



lighter shading (these flows are also in red in the online ver- 
sion). For illustrative purposes, we only show the flows for 
when the temperature is not sufficiently hot for photodissoci- 
ation of '^N. Finally, by convention and to avoid cluttering 
the plot, we do not show flows into or out of p, n, and '*He. 
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Fig. 2. — Flows during a constant-pressure self-heating burn at p = 
lO' gcm"^ with an initial composition {X12, Xie, X22I = {0.3,0.7,0.0). We 
integrate over the period of the bum with T < 6 X 10** K. Each row contains 
the isotopes of a particular element Z, with the columns containing different 
neutron numbers N. The width of the arrow is proportional to the magnitude 
of the flow, with only those flows having magnitude > 0.01 ■ max(F) being 
shown. The primary '^C + reactive flows are indicated with a lighter 
shading and are in red in the online version. 
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Fig. 3. — Same as in Fig. [2] but at p 
^'Ne is now tl 
(p, Q"), and (n, y) reactions. 
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The reaction 



is now the dominant destroyer of -^Na instead of (p,y). 



To facilitate comparisons between different runs, we define 
a normalized flow as 



F(i ^ ;■) = 



F{i ^ j) 



FC'-C ^ 23Na) + FC^C ^ 20Ne)' 



(4) 



At a density of 10^ gcm"^^ we find ^(i^C ^ ^3^^^) ^ Q 43 
andf'('^C ^°Ne) = 0.57 , which reflects the branching ratio 
(ICa ughlan & Fowle^l 19881) . In agreement with the arguments 
in §|2l most of the protons liberated by '-C('^C, /:')^''Na cap- 
ture onto i^c, with FC'C '^N) ^ FC^N '^Q = 0.20. 
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Note that if all p were to capture onto '^C, we would have 
/■(12c ^ ^^N)/FC^C ^^Na) = 0.5, because the re- 
action '^C('^C, /7)^^Na consumes two '^C. The o'-particle 
released by the '^C('^C, a)^'^Ne reaction captures onto '■'C 
to form ""O and a neutron, which in turn destroys another 
'^C via '2C(«,7)i^C. The flow FC^C -» '^Q = 0.26 
and nearly matches the number of a-particles produced by 
'^C('^C,Q')^''Ne; «-captures onto ^■'Na and ^°Ne account for 
the difference. At p - 10^ gcm"^ the only electron cap- 
tures are onto '■'N, so AY^ - -FC^N '■'C). Dividing by 
AFi2 = -[FC^C ^^Na) + FC^C ^°Ne) + PC^C 
i^N) + FC^C "C)] gives dYe/dYn = 0.14 (Table[I]i. This 
is slightly less than the estimate of 1/6 (§ I2.2| | because of the 
lower branching ratio of '^C — » ^^Na. 

At p = 3 X lO'' gem""*, "■'Na is consumed by the reaction 
^■'Na(e", Ve)--'Ne rather than by p- or «-capture (cf. Figs.[2]and 
[3]l. The flow is much larger than that from '^N, because the 
reaction ^^Ne{p,n)-^Na competes for p and produces more 
23Na. Indeed, ^('^C ^ '^''C) = 2.5FC^C "N) because 
of extra « produced by the reaction -^Ne{p, «)^-'Na. There 
is an additional contribution from electron captures onto '^F 
produced via '^0(/7,y)'^F, but this flow is only about 4% of 
the -^Na(e' , Ve)-^Ne flow. At both densities, ^^Ne plays a 
small role in reducing Y^ which we verified by comparing 
the flows for a burn with X22 = 0.06 with those for a flow 
with X22 - 0.0 (Fig. |3]l. We find that, for a burn starting at 
p = 3.0x10'' gcm-^,F(^^Ne ^^Na)/FC^C '^N) = 0.08 
and PC^Ne ^^Ne)/PC^C ^^C) = 0.05. Although 
/lpNe(/7,7)23Na] < A[^^C(p,yy^N] at T < 6 x 10** K 
(Fig. [U, the abundances are in ratio Y22/Y12 = 0.11. Be- 
low the electron capture threshold for ^^Na, the reaction 
^^Ne{p, y)^-'Na will cause a slight reduction in dYddYn equal 
to the ratio of ^(^^Ne ^ ^^Na) to ^('^C ^ ^^N). 

Finally, we confirmed the lack of an s-process flow (§ l2.3l l 
by performing a run with nuclides from Ne to ^^ Fe present 
at up to 3 times their solar abundances (lAnders & Grevess^ 
Il989l) . The threshold for electron capture onto ^^Fe is 1.5 x 
10^ gcm"^, and so at higher densities carbon ignition occurs 
in a more neutron-rich environment. We therefore start the 
calculation by artificially suppressing the strong interactions 
and allowing the mixture to come into y6-equilibrium. We then 
turn on the strong reactions and let the runaway commence. In 
all cases the heavier nuclides did not have a substantial impact 
on the reactive flows. The value for dYf,/dYi2 is somewhat 
larger than 0.3 at densities p > 1.2 x 10^ gcm"^, the thresh- 
old density for ^^Mg(e", Vf,)^^Na(p, «), because of the reac- 
tions ^'*Mg{p,y)^^Al{e',Vef^Mg(e',Ve)^^Na. Of these two 
captures, Al is y6^-unstable, and hence the electron capture 
is fast enough to proceed throughout simmering; the capture 
onto ^^Mg has a timescale, at p = 2 x lO'^ gem"-', of » 900 s 
and will therefore freezeout during simmering, just as cap- 
tures onto ^^'Na freezeout. 

Our runs span a range of initial densities, from lO** g cm""* 
(for which the electron Fermi energy is too low to induce elec- 
tron captures onto ^^Na) to 6 x 10'^ gcm"^, which represents 
an extreme c ase for accretion on to a cold, initially massive 
white dwarf (iLesaffre et al.l 120061) . In all cases we took the 
initial '^C mass fraction to be X12 - 0.3. As noted above, 
the mass fraction of ^^Ne would have to exceed that of '^C 
to change the nucleosynthesis during simmering appreciably. 
At higher densities, p-captures onto -^Mg can also play a role 
(§ 13. It . but our results will not change appreciably so long as 
X12 is not substantially less than 0.3. 

The calculation, being a single reaction network integra- 



TABLE 1 

Change in electron abundance per carbon consumed 
during the pre-explosion convective burning 



composition'^ 


density 
10' gcm-3 


aigiar 12 


Li "'i/"'l2 


X22 = 0.00 


1.0 


0.136 


0.340 




3.0 


0.297 


0.340 




6.0 


0.302 


0.342 


X22 = 0.06 


1.0 


0.125 


0.347 




3.0 


0.301 


0.344 




6.0 


0.305 


0.346 




1.0 


0.130 


0.361 




3.0 


0.349 


0.370 




6.0 


0.355 


0.380 



" In all cases the initial mass fraction of '^C is 0.3. 



tion, does not incorporate the effects of mixing in the white 
dwarf core. Our focus here is to elucidate the reactions that 
set Yg. These calculations do not determine the total amount 
of carbon consumed (although see § l3.3l for an estimate) or the 
total mass of processed material that lies at a density greater 
than the electron capture threshold. We list our one-zone re- 
sults in terms of the change in electron abundance per carbon 
consumed, dY^ldYn. Table[T]summarizes our numerical find- 
ings of dYJdYi2 for densities 10^, 3 x lO'^, and 6 x 10^ g cm"-', 
for ^^Ne mass fractions X22 - and 0.06, and finally a run 
(denoted as 3Zq in Table [1]) with elements heavier than ^"^Ne 
present at 3 times solar abundance. We use this value of 3Zq 
as representing a rough upper limit based on the ^ 0.5 dex 
scatte r in [Fe/H] present in local field stars (iFeltzing et al.l 
I2OOI . 

3.2. The effective q-value of the '^C + '^C reaction 

The scope of this work is to elucidate the nuclear reactions 
that occur during the pre-explosion simmering, and including 
their effects in large-scale hydrodynamics simulations of the 
entire white dwarf is advisable. As an aid to such simula- 
tions, we present fits for the carbon depletion rate and effec- 
tive heat deposition, which improve on previous expressions 
dWooslev et al.ll2004 . Since the reaction chain is controlled 
by the reaction '^C-h'^C, we write the rate of carbon consump- 
tion, dY\2ldt as being proportional to the thermally averaged 
cross-section, A - Na{o-v} 

This definition is such that the quantity in parenthesis is the re- 
action rate per pair of '^C nuclei and M12 = 2 if the only '^C- 
destroying reaction present is '^C + '^C. One can determine 
M12 from summing the normalized reaction flows (eq. Q) out 
of 12c. 

In a similar fashion, we can define the effective heat 
evolved, q^-ff, per reaction '^C + '~C by the equation 

e^^eifA^A^^J'nP^)- (6) 

With this definition, one has e = iqesNp^IMn) x dYn/dt; 
this may be compared with equation (|2]i to relate ^eff to the 
binding energy of the nuclei. To compute these quantities, 
we integrated the reaction network over a grid of p and T, 
both of which were held fixed for each run. We found in all 
cases that the instantaneous energy generation rate s would, 
after some initial transient fluctuations, settle onto a constant 
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value until a significant (> 10%) depletion in '^C had oc- 
curred. We used this plateau in e to obtain ^eff and M12. For 
the densities of interest, the values of both ^eff and M12 ob- 
tained this way are nearly independent of temperature. We 
find M12 - 2.93; this value is accurate to within 2% over all 
our runs. The value of (^eft increases slightly with density, but 
is nearly constant over the temperature range of interest. We 
find that q^s = 8.91 MeV, 9.11 MeV, and 9.43 MeV for 
p = 10^ gem"-', 2 X 10^ gem"-', and 3 x 10'' gem"-', respec- 
tively. At each of these densities, the quoted value of q^ff is 
accurate to within 3% over the temperature range of 3 x 10*^ K 
to 7 X 10*^ K. 

For use in large-scale hydrodynamic models, one first com- 
putes the screened the rmally averaged cros s-section A (for the 
most recent rate, see iGasques et al.ll2005h at the thermody- 
namic conditions of a given cell. Combining A and the cell's 
carbon abundance Y12 with our estimates of M12 and g'eff, one 
computes dYnldt and e from equations (|5]) and In ef- 
fect, this procedure incorporates the results of a large reaction 
network and careful treatment of the detailed nuclear physics 
into simple expressions. We caution that these fits were ob- 
tained in the regime -AF^ < 0.003, for which proton cap- 
tures onto ^^Na are not competitive with proton captures onto 
'^C. Note that the heat released per '^C nucleus consumed 
is qeff/Mi2 ~ 3.1 MeV, slightly higher than our simple es- 
timate (§ [2li, and also somewhat higher than those used by 
iPiro & BildstenI (l2008h . This is because of our inclusion of 
heating from the '■'N(e", Vp)'-'C and -■'Na(e", V(.)^-'Ne reac- 
tions. Our estimate of the heat evolved is less than that com- 
puted under the assumption that the pr oducts of '^C burnin g 
are a 3:1 -°Ne:^**Mg mixture (see, e.g. IWooslev et al.ll2004b . 
which releases 5.0 MeV per '^C nucleus consumed*. Our net 
heating rate, per '^C-h'^C rea ction, is about 10% less than that 
used by lWoosley et al.l(l2004h . but we effectively consume w 3 
'^C nuclei per reaction. 

The change in electron abundance is related to A via 



dt 



dV, 
'dYT-, 



1 . 
-Y 

2 ■ 



pA 



(7) 



where dYe I dY 12 is taken from Table [T] Finally, we may com- 
pute the change in the mean atomic mass, {A), as a function of 
carbon consumed. On differentiating (A) = {YjidYildY\2y^ , 
and substituting equation (|5] we have 



d{A) 
dt 



{AfMn\^Yl^A 



V ^ 
4^ dYn 



In previous sections, we evaluated the heating and neutron- 
ization of the white dwarf core in terms of the rate of '^C 
consumption. We now estimate the net amount of '^C that is 
consumed in raising the white dwarf temperature and eval- 
uate the temperature at which electron captures onto ^■'Na 
freeze out. In the one-zone isobaric calculations, using equa- 
tions © and I©, we have dT/dYu - {qesNAjliMnCp), so 
that Ar * 0.15 GK(AXi2/0.01), where AX12 is the change in 
mass fraction of '^C. The change of '^C abundance required 
to raise the temperature from 3 x 10^ K to 8 x 10^ K is then 
lAFnl « 2.8 X 10"l This is about 11% of the available '^c, 
for an initial '^C mass fraction X\2 - 0.3. Figure |4] shows 
the decrement of the e" abundance, (f = 0) - Yi,{t) - -AY^ 
as a function of '^C consumed, Yi2{t = 0) - Fi2(0 - -Ayi2. 
Note that we are plotting the decrement in abundance. The 
calculation was started at an initial density and temperature 
p = 3.0 X 10'' gcm"^ and r = 1.9 X 10^ K. As tu shortens, 
the electron captures onto -^Na "freeze out" and dYg/dY^ de- 
crease to X! 1/6. When -AF12 > 0.003, the abundances of 
^^Na and -•'Ne have increased sufficiently that they compete 
with to consume pr otons, and thereby halt the neutroniza- 
tion, in agreement with lPiro & BildstenI d2008l) . 



o 
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12 3 4 

[yi2(/ = 0)-Fl2(0]xl03 

Fig. 4. — Change in electron abundance, ^^(f = 0) - Yg{(), as a function of 
carbon consumed, }'i2(? = 0) - Fi2(f)- The break in the slope, at Fi2(0) - 
Y\2 « 2 X 10"^, occurs when the heating timescale ?h = Cp/e becomes 
less than the timescale for electron capture onto ^^Na, which is » 2700 s at 
p = 3.0 X lO' gcm~^. We indicate this point with the thin vertical line. To 
guide the eye, the short dotted line indicates a slope of \^ . 



where the quantity YjidYijdY\2 is computed from the flows 
(eq. HI) and is listed in Tabl e [U From the sirnple d escription 
of the reactions (§|2l see also lPiro & Bildstenll2008h we have, 
for every six '^C destroyed, one each of '■'C, '*0, ^°Ne, and 
23Na (or 23Ne), so that Y^id^ild^n ~ f This results in a 
smaller change in (A) than would result from burning '^C -H 
'^C to a 3:1 ^°Ne:^^^Na mixture (in that case, YjidYi/dYn - 
0.42). We advocate using equations (I5]l-(l8]l and the computed 
values (Table [U of dY^ldY^ and 2, c/i',/c/i'i2 in numerical 
simulations of simmering. 

3.3. Heating of the white dwarf core and the end of 
simmering 



' Note that in eq. (1) of lWooslev et alj <2004) . the factor of 
into their quantity /Ii2,i2- 



is subsumed 



As noted bv lPiro & BildstenI (|2008|) . a one-zone model will 
overestimate the heat required to raise the central tempera- 
ture by a given amount, and hence overestimate the amount 
of '^C that must be consumed during simmering. We per- 
form a calculation similiar to theirs: starting with an isother- 
mal white dwarf with a given central density and with a tem- 
perature set by equating heating from the '^C + '~C reaction 
with neutrino losses, we then raise the central temperature to 
8 X 10** K, keeping the total white dwarf mass fixed, and fol- 
low an adiabatic temperature gradient until we intersect the 
original isotherm at radius rconv, which we then follow to the 
stellar surface. We compute the total stellar energy, gravita- 
tional and thermal, in both cases, and take the difference to 
find the heat required to raise the central temperature of the 
star to 8 X 10** K. The temperature of 8 x lO** K is chosen 
as a fiducial temperature representing the point at which the 
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heating of a fluid element proce eds faster than the convective 
timescale (IWoosley et al.ll2004l) . The evolution of the white 
dwarf core is not exactly isobaric: the expanding convective 
zone heats the white dwarf. As the entropy of the white dwarf 
increases, it expands and reduces the core pressure. For an 
initial central density pinit = 3.0 x 10^ g cm"^, we find that in 
raising the central temperature from 3 x 10** K to 8 x 10^ K 
the radius expands by a factor of 1 . 1 and the central pressure 
decreases to 0.59 of its initial value. 
For initial central densities pinit = 10"^ gem"-', 3 x 

,-3 



10^ gcm^^*, and 6 x lO** gcm'\ corresponding to white 
dwarf masses of 1.35 Mq, 1.38 Mq, and 1.39 Mq, the ini- 
tial temperatures defined by the onset of thermal instability 
are 3.8 x 10** K, 1.9 x 10** K, and 1.0 x 10** K, respec- 
tively. The energy required to raise the central temperature 
to 8 X lO** K is E, = 2.11, 4.12, and 5.58 keV per nu- 
cleon, respectively. When the central temperature has reached 
8 X 10^ K, the masses of the convective zone for these three 
cases are 0.69 Mq, 1.10 Mq, and 1.29 Mq, respectively. 
The spatial extent of the convective zone, rconv ~ 1000 km 
for the three cases, is in agreement with the findings of 
iKuhlen et al.l (l2006h . We checked our computation of Ec by 
computing E'^ = LCp [Tf,„^\ - Tinitiai] dM, as was done by 
iPiro & BildstenI (l2008h . Both methods give comparable es- 
timates, but slightly underestimates the change in energy 
(by »J 10%), because it does not account for the expansion 
of the white dwarf. Neglecting the change in 7^ as the white 
dwarf heats introduces a small correction to E^-: for an adi- 
abatic 1.38 Mq white dwarf with a central temperature of 
8 X 10'^ K, a reduction in by 1.66 x 10"^^ reduces E,- by 
only 3.3%. 

If the white dwarf were entirely mixed, raising the cen- 
tral temperature 8x10^ K would require, for the three pinit 
cases here, that lAFnl = 1.36 x 10"^ 1.66 x 10"^ and 
1.87 X lO"-', respectively. Because the changes in composi- 
tion are only mixed over the convective zone, the decrement 
in Yi2, and hence Yf,, is more pronounced there. Using our 
estimate of rconv, we estimate that within the convective zone 
lAFnl = 3.0 X 10"\ 2.1 x 10"^ and 2.01 x 10"^ respectively, 
forpinit/ClO"^ gcm-3) = 1.0, 3.0, and 6.0. Should the radial 
extent of the convective zone be smaller tha n our estimate, for 
example because of convective Urea losses ( Stein et aO 19991: 
IStein & Wheeled 120061: iLesafire et all 120051) . the abundance 
of '^C will be further reduced in the white dwarf core. A 
lower '^C abundance reduces the laminar speed of the flame 
launched at the end of simmering (iTimmes & Woosley|[T992l: 
IChamulak et aDl2007h . 

Finally, we estimate at what point the convective mixing 
timescale becomes shorter than the timescale for electron cap- 
tures onto ^''Na. Using our adiabatic temperature-gradient 
white dwarf models, we compute the typical convective ve- 
locity v,n,s, and hence a characteristic turnover time fconv - 
rcnnv/vrms, using the mixing le ngth formalism (see the dis- 
cussion in IWooslev et al.ll20()4l) with the total luminosity and 
evaluating thermodynamical quantities at their central values. 
Figure |5] shows fconv {thick lines) for the three cases of pinit 
considered above, as well as the electron capture timescale, 
fee for ^■'Na [thin lines) for those cases with a density above 
threshold. To estimate the effect of the ^^Na/^^Ne pairs on the 
convective zone, we evaluated the mass fraction of these pairs 
for the case pinit - 3.0 x 10'' gcm^^* when the temperature 
had risen to T - 3.5 x 10** K and fee ~ feonv (Fig- E] solid 
line). At this point the convective core has a mass 0.5 Mq 



which is comparable to their calculation. We estimate, from 
the energy required to heat the white dwarf to this point, that 
the mass fraction of ^^'Na/^^'Ne pairs in the convective zone 
will be X23 = 0.004 at this time. We note that such a large 
number of Urea pairs will ha ve a dramatic effect o n the prop- 
erties of the convection zone (iLesaffre et al.ll2005h . As IS evi- 
dent from Figure|5] there is a range of temperatures for which 
fconv < fee < ^H- In this region, the effective dYeldY\2 will de- 
pend on the fraction of mass with densities above the capture 
threshold and on the effects of Urea losses on the convective 
flows, but will still be larger than the minimum value set by 
'"'N(e", V(,)''C. Incorporating the effects of the neutrino lu- 
minosity from the ^-^ Na/^-^Ne reactions in this re gime is nu- 
merically challenging (ILesaffre et al.ll2003ll2006l) and beyond 
the scope of this work; for now, we just note that this is the 
primary uncertainty in determining the amount of con- 
sumed during the pre-explosive phase, and a better treatment 
is needed. 
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Fig. 5. — Mixing lengtli convective dmescales (thick lines) and electron 
capture timescales (for densities above tliresliold) for ^^Na (thin lines), as 
a function of the central temperature during the simmering. Three initial 
densities are shown: lO' gcm"^ (dotted line), 3 X lO' gcm"^ (solid lines), 
and 6 X lO' gcm"-^ (dashed lines). The electron capture timescales increase 
as the white dwarf heats because the density decreases during simmering. 



4. DISCUSSION AND CONCLUSIONS 

Using a nuclear reaction network coupled to the equation 
for self-heating at constant pressure (eq. ID-lEl), we have in- 
vestigated the change in Ye induced by electron captures on 
nuclei produced by '^C fusion during the pre-explosion sim- 
mering of the white dwarf . We confirm that ther e is a maxi- 
mum Ye at flame ignition dPiro & BildstenI l2008h . We quan- 
tified the role of ^^Ne and other trace nuclides in setting the 
change in electron abundance with '^C consumption by us- 
ing a full reaction network, and we included the heating from 
electron captures into an excited state of '^C. We gave simple 
formulae (eq. ||5|-|l8l) for the energy generation rate, the rate 
of change in electron abundance, and the rate of change in the 
mean atomic mass to include the detailed nuclear physics into 
large-scale hydrodynamical simulations. 

Our estimates of the maximum Ye at the time of t he exp lo- 
sion are roughly similar to those of lPiro & BildstenI ( 120081) . If 
we neglect the effect of thermal neutrino losses on the evolu- 
tion of the white dwarf, then Ye is reduced by 2.7 x 10 '*-6.3 x 
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10 within the convective zone. This reduction in Yg depends 
predominantly on the amount of '^C consumed prior to igni- 
tion. The electron captures during simmering reduce Yg be- 
low the value set by neutron-rich -^Ne inherited from core He 
burning by the white dwarf's progenitor star. Reducing Y^ 
in the explosion depresses the yield of ^ ^Ni and increases the 
amount of ^^Fe and ^^Ni synthesized dlwamoto et al.l IT999t 
iTimmes et al.ll2003l) . even in the absence of further electron 
captures onto the Fe-peak isotopes in nuclear statistical equi- 
librium (NSE) in the densest portion of the white dwarf. As 
a result, any correlation between host system metallicity and 
white dwarf peak luminosity will be weakened for Z < Zq (for 
which the reduction in 7^ due to captures during simmering is 
greater than the change due to initial white dwarf composi- 
tion). 

To illustrate how the simmering electron captures affect 
the light curve , we r econstruct the comparison made by 
[Gallagher et al.l ( l2005h . who compiled a set of SNe la with 
measured AmisiB), defined as the change in B over 15 days 
post-peak, and host galaxies with measured abundances of 
oxygen to hydrogen, denoted O/H. We construct an expres- 
sion for M56, the mass of ""^Ni produced in the explosion, that 
depends on dYddYn (Table [U, AY\2, and the host galaxy 
composition (we assume that the white dwarf has the same 
O/H ratio as the galaxy). The trace nuclide that predomi- 
nantly sets Ye in the white dwarf is ^^Ne, which traces the 
aboriginal abundance of CNO nuclei in the main-sequence 
star from which the white dwa rf evolved. We ther efore re- 
cast the linear formula for M56 (iTimmes et al.ll200"3h in terms 
of O/H. For simplicity, we fix the 'H:^He ratio, as well as the 
ratio of heaver elem ents to '^O, to their solar system values 
dAsplund et al.ll20()5l) a nd neglect correction s from the change 
in [O/Fe] with [Fe/H I ([Ramfrez et al.ll200'7l) and the increase 
in [N/O] with [O/H] (ILiang et al.ll2006l) . With these assump- 
tions, the mass of ^^Ni ejected in the explosion is 



2.0 - 



M56 = M56,0 



,0\ dYe 



(9) 



where Mse.o is the total mass of NSE material synthesized 
at densities where electron captures during the explosion are 
negligible, dY^dYn is the change in electron abundance with 
carbon consumption given in Table [1] and -AKn is the total 
amount of '^C consumed during the pre-explosion convective 
phase. 

Relating this to an observable such as Ami5(B) requires an 
explosion model that follows the radiation transfer. Our inter- 
est here is to isolate how Mse changes with the CNO abun- 
dance of the progenitor when the other parameters, such as 
the ejecta kinetic energy and the total mass of iron-peak nu- 
clei, are held fixed. This is important, because the relation 
between light curve width and pe ak brightness depen ds on 
these other parameters as well (see lWooslev e t al. 200 7t for a 
lucid discussion). We use the model M070103 dWooslev et al.l 
l2007h . in which the total mass of iron-peak ejecta is 0.8 M©: 
of that, the innermost 0.1 Mq is stable iron formed in situ 
from electron captures during the explosion, with the remain- 
der being a mix of radioactive ^^Ni and stable Fe. Note 
that this model follows the peak lu minosity-light curve width 
relatio n (iPhillipsI Il993h . whereas iMazzali & Podsiadlowskil 
(I2OO6I) suggest that varying the ratio of ^^Ni to stable Fe may 
create dispersion abou t this relation. 

Using the model of ' Woosley et al.l (l2007l) . we set M56.0 = 
0.7 Mq, use equation (|9]l to compute M5 6 as a function of O/H 
for different l^Y\2, and interpolate from I Woosley et all ( l2007i 
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Fig. 6. — Conelation between 12 + log(0/H) and Am 15 (B) induced by 
the ele ctron abundance. The points in the plot are taken t'rom lGallagher et aU 
120051) . We show the 'maximal' case for the reduction in Y,, during simmering 
{dashed line), in which dYgldY\i = 0.30 in the early stages of the bum, 
which is an upper limit to the derivative. In this case AFf = -7.5 X 10"^. For 
comparison, we also show the case in which there is no reduction of Y,, during 
simmering (solid line). Not e that the relation betw een M56 and Ami5(B) is 
shallower than that used bv lGallagher et al] 120051 dotted line); see the text 
for an explanation. To the left of the vertical bar, the decrease in Msg due 
to electron captures during simmering will exceed that due to enrichment by 
^^Ne, and hence any correlation between O/H and Ami5(B) will be masked. 

Fig. 22) to find ^nii(B). Figure |6] displays this result. We 
plot here a maximal case {dashed line) with dY^/dYn = 0.30 
in the initial part of the simmer, appropriate for the one-zone 
calculation with electron captures onto ^^Na (Fig.|4]i; for this 
case AYe = -7.5 x 10""*. For comparison, we also plot a 
case {solid line) with AY^ - during simmering. This gives 
a sense of how large the variation in Yg might be. The ver- 
tical bar indicates the value of 12 + logj,)(0/H) at which 
the change in Amis{B) from the ^^Ne abundance equals that 
from the electron captures during simmering for this maximal 
case. To the left of this curve the linear correlation between 
^-Ne abundance and Amis{B) will be masked by variations 
in the simmering of the white dwarf. For comparison, we 
also sh ow the data from the c ompilation of [Gallagher et alj 
(l2005h and lHamuv et al.1 (|2000|) and p lot the relation b e tween 
12 -I- log(0/H) and Amii{B) used bv iGallagher et all (I2005L 
dotted line). This trend is much steeper than our finding. The 
difference is due to how the ^^Ni mass was varied; whereas 
the models we use ( IWooslev et alj|2007l) hol d the kinetic en- 
ergy a nd total mass of iron-peak ejecta fixed. [Gallagher et alJ 
(l2005h based their peak b rightness on delayed detonation 
models ( iHoflich et alj2002h for which a variation in ^^Ni also 
produced a change in the relative amounts of iron-peak and 
intermediate mass-elements, as well as a different explosion 
kinetic energy. 

It is evident from Fig. |6] that the scatter in the data points 
is larger than the expected trend due to progenitor compo- 
sition, especially at sub-solar metallicities. Both AF^ and 
dYeldYi2 depend on the central density, which is not ob- 
viously correlated with metallicity, and hence the correla- 
tion between peak brightness and O/H will be masked by 
differences in the pre-explosion simmering. Indeed, if the 
variation in AYe were as large as the two cases we plot in 
Fig.|6] then variations in ^*'Ni would be determined more by 
AY\2 than by stellar composition for galaxies with sub-solar 
O/H. There is a general trend that SNe la are systematically 
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brighter in galaxies with active star-formation ( 


Hamuv et al. 


200C 


;lGallagher et al.l2005HSullivan et al.l200(l 


Howell et al. 


2007 


). ISuUivan et al.1 (|2006|) showed that the SNe la rate 



increases linearly with the specific star formation rate, and 
that SNe la associated with actively star forming galaxies 
were intrinsically brighter than those associated with pas- 
sive galaxies. Although many of these passi ve galaxies are 
more massive, and hence more metal-rich dTremonti et alj 
120041) . the observed scatter in SNe la peak brightnesses re- 
mains much larger tha n the expected trend with metallicity 
jPiro & BildstenlboOSi : D. A. Howell 2007, private communi- 
cation). This suggests that the correlation with the chemical 
abundances of the host galaxy is a secondary effect in setting 
the peak brightness of SNe la. 
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